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^ '■ Abstract 

j^ [ The dynamical nucleus-nucleus potentials for some fusion reactions are investigated by using 

the improved quantum molecular dynamics (ImQMD) model with different sets of parameters in 

which the corresponding incompressibility coefficient of nuclear matter is different. Two new sets 

of parameters SKP* and IQ3 for the ImQMD model are proposed with the incompressibility coef- 

Q ' ficient of 195 and 225 MeV, respectively. The measured fusion excitation function for ^^0+^'^^Pb 

P3 ' and the charge distribution of fragments for Ca+Ca and Au+Au in multi-fragmentation process 

can be reasonably well reproduced. Simultaneously, the influence of the nuclear matter incom- 

vQ I pressibility and the range of nucleon-nucleon interaction on the nucleus-nucleus dynamic potential 

(^ ' is investigated, 
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I. INTRODUCTION 

The synthesis of super-heavy nuclei and heavy- ion fusion at deep sub-barrier energies have 
attracted a great deal of attention in recent years [1-8]. The calculation of nucleus-nucleus 
potential especially at short distances is of crucial importance for these studies. Several 
static and dynamic models have been proposed for calculating the nucleus-nucleus potential 
[9-15]. The static potentials are usually obtained by some empirical formulas or based 
on the double-folding concept and the sudden approximation. To consider the dynamical 
process in fusion reaction microscopically, some microscopic dynamics models, such as the 
time-dependent Hartree-Fock (TDHF) model [13, 14] and the improved quantum molecular 
dynamic (ImQMD) model [16, 17] have been developed. The ImQMD model is a semi- 
classical microscopic dynamics model and is successfully applied on intermediate-energy 
heavy-ion collisions and heavy-ion reactions at energies around the Coulomb barrier [15, 
17, 18]. In Ref. [15] the extended Thomas- Fermi approximation is adopted for calculation of 
the dynamical nucleus-nucleus potential at short distances based on the obtained dynamical 
densities of the reaction system from the ImQMD simulations. The energy dependence 
of the dynamical nucleus-nucleus potential was observed. In addition to the influence of 
the incident energy on the nucleus-nucleus potential due to the dynamical evolution of the 
densities, the influence of the incompressibility coefficient of nuclear matter on the nucleus- 
nucleus potential and the fusion cross sections are also investigated in Refs. [19, 20] according 
to the double-folding calculation with the M3Y interactions. 

To investigate the influence of nuclear equation of state on the djTiamical nucleus-nucleus 
potential, we will study the nucleus-nucleus potential at short distances with the ImQMD 
model by adopting different sets of parameters. The corresponding nuclear equation of state 
for these different sets of parameters is different, with which we attempt to understand 
the influence of the incompressibility coefficient of nuclear matter on the dynamical fusion 
potential. The structure of this paper is as follows: In sec. II, the ImQMD model is 
briefly introduced. In sec. Ill, two sets of parameters IQ3 and SKP* are proposed for 
the ImQMD calculation, and the fusion reactions ^^O+^^^Pb and ^^Ca+^°*Pb at energies 
around the Coulomb barrier and the reactions Ca+Ca and Au+Au at incident energy of 
35MeV/nucleon are also studied for testing IQ3 and SKP*. In addition, the dynamical 
nucleus- nucleus potential inside the Coulomb barrier is simultaneously investigated. Finally 



the conclusion is given in Sec. IV. 

II. THE IMPROVED QUANTUM MOLECULAR DYNAMICS MODEL 

In the ImQMD model, as in the original QMD model [21], each nucleon is represented by 
a coherent state of a Gaussian wave packet. The density distribution function p of a system 
reads 

^_ 1 r (r---r.V\ 

(1) 
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where ar represents the spatial spread of the wave packet. The propagation of nucleons is 

governed by Hamiltonian equations of motion under the self-consistently generated mean 

field, 

dH . dH .. 

where r, and Pi are the center of the i-th wave packet in the coordinate and momentum 

2 

space, respectively. The Hamiltonian H consists of the kinetic energy 7" = ^ ^ and the 

i 

effective interaction potential energy U: 

H = T + U. (3) 

The effective interaction potential energy is written as the sum of the nuclear interaction 
potential energy U\oc = J V\oc{r)dr and the Coulomb interaction potential energy f/coui which 
includes the contribution of the direct and exchange terms, 

U = Uloc + f/coui. (4) 

Where Vioc(?") is the potential energy density that is obtained from the effective Skyrme 
interaction and taken to be the same as that in Ref. [17]: 

Vice = f ^ + ^^ + f^(Vp)^ + g/-^ + ^[P'- UVpr]6^ (5) 

2 po 7 + 1 Po 2po Po' 2po 

where S = {pn — Pp)/{Pn + Pp) is the isospin asymmetry. To describe the fermionic nature 
of the N-body system and to improve the stability of an individual nucleus, the phase- 
space occupation constraint method [22] and the system-size-dependent wave-packet width 
(Jr = ctq + o'lA^/'^ fm [16] are adopted. The parameter sets adopted in this work are shown 
in Table I. 



TABLE I: Model parameters adopted in this work. 



Parameter a /3 7 (7sur 9t V Cg Hg po ctq cji 

(MeV) (MeV) (MeVfm^) (MeV) (MeV) (fm^) (fm-^) (fm) (fm) 

IQ2 -356 303 7/6 7.0 12.5 2/3 32.0 0.08 0.165 0.88 0.09 

SKP* -356 303 7/6 19.5 13.0 2/3 35.0 0.65 0.162 0.94 0.018 

IQ3 -207 138 7/6 18.0 14.0 5/3 32.0 0.08 0.165 0.94 0.018 

III. RESULTS 

In this section we first briefly introduce the parameter sets IQ3 and SKP*. Then we test 
the new parameter sets through fusion reactions and heavy-ion coUisions at intermediate 
energy. Finally we investigate the dynamical nucleus-nucleus potential inside of the Coulomb 
barrier. 

A. New parameter sets SKP* and IQ3 

According to the properties of nuclei at ground state and the knowledge of nuclear in- 
compressibility, a number of Skyrme forces, such as SkM* [23], SKP [24] and SLy4 [25] were 
proposed in recent decades. With the proposed parameter sets, the Skyrme energy-density 
functionals have been successfully applied on the studies of nuclear structure, fusion reaction 
and neutron star, etc. Based on the parameters of Skyrme forces, the ImQMD parameters 
can be directly obtained as those done in Ref. [18]. Considering the range of nucleon-nucleon 
interaction in the ImQMD model which is represented by a gaussian wave-packet, the pa- 
rameters of the ImQMD model in this work are re-adjusted for studying the fusion reactions 
in which the stability of an individual nucleus plays a role for a reliable description of the 
dynamical process. 

The knowledge of nuclear equation of state at densities around the normal density po 
is helpful to constrain the model parameters. In the ImQMD model, the incompressibility 
coefficient K^o of symmetric nuclear matter at po is expressed as 



Koo = 9p| 
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with Cfc = ^l(^)^''^ = 75.0 MeV fm^. ^ = co/cfc is a correction factor for the kinetic energy 
of a nuclear system when applying the extended Thomas- Fermi (ETF) approximation which 



is roughly expressed as 

Er - CO 5: pT + v^ e f^P^^ + ^-^ (7) 

in this model. The expressions of pi and pij are given in Ref. [15]. The coefficients cq, ci 
and C2 can be determined by the kinetic energies T of nuclei at their ground state. ^ = 1 is 
for the idealistic Fermi-gas. Considering the range of realistic nucleon-nucleon interactions 
which is described by gaussian wave-packets in the ImQMD model, we approximately set 
C, = 0.4 ~ 0.6 for a reasonable description of the properties of nuclei at ground state and the 
stability of an individual nucleus. The corresponding energy per particle E/A at po reads 



E{p) 



A 



= ^CkpT + ^ + ^+9r. (8) 

P=Po 2 7 + 1 



Although experimental and theoretical investigations on the nuclear equation of state suggest 
that KoQ ~ 230 MeV, E/A ^ —16 MeV around the saturation density which is ~ 0.16 
fm~^ [26-28], the uncertainty of nuclear equation of state still causes some difficulties for a 
unambiguous determination of the model parameters. 

To investigate the influence of the nuclear matter incompressibility coefficient on fusion 
reactions, we attempt to propose two new sets of parameters, SKP* and IQ3 for the ImQMD 
simulations. The parameter set SKP* is generally determined based on the Skyrme force 
SKP [24] in which the parameters Qr, Cq and ai are adjusted for an appropriate description 
of nuclear properties at ground state and the fusion reactions. For both 1Q2 and SKP*, the 
corresponding values of Koo are the same but the wave-packet widths ar are different, which 
is useful for exploring the influence of the interaction range of nucleons and the flnite-size 
effect of nuclei. To explore the influence of the nuclear matter incompressibility on the 
fusion reactions, we also construct the parameter set IQ3 in which the wave-packet widths 
ar are the same as those in SKP* but the incompressibility coefficient K^o = 225 MeV is 
obviously larger than that in SKP* and IQ2 (see Table II). The corresponding kinetic energy 
coefficients for different sets of parameters are listed in Table II. 

B. Tests for IQ3 and SKP* 

With the parameter sets IQ3 and SKP*, the time evolutions of the binding energies 
and nuclear radii for a number of nuclei have been checked. We flnd that an individual 



TABLE II: Kinetic energy coefficients in the ETF approximation and incompressibility coefficient 
for parameter sets IQ2, SKP* and IQ3. 



Parameter 


Co (MeV fm^) 


Cl 


(MeV /m^) 


C2 


(MeV) 


Koc (MeV) 


IQ2 


41.2 




4.8 




-1.0 


195 


SKP* 


38.9 




0.1 







195 


IQ3 


43.3 




0.34 







225 
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FIG. 1: (Color online) Time evolution of the density distribution for the fusion reaction Ca+^^^Pb 
at an incident energy of -Ec.m. = 200 MeV with IQ3. 

nucleus can remain stable for several thousands fm/c without spurious nucleon emission. 
Simultaneously, some fusion reactions are investigated for testing the parameter sets. Fig. 1 
shows the time evolution of the densities for the fusion reaction ^^Ca+^°^Pb at an incident 
energy of -Ec.m. = 200 MeV with the parameter set IQ3. From the time evolution, one sees 
that the central densities of nuclei in the reactions are reasonable. In addition, we note that 
the surface diffuseness of nuclei at neck region increases when the neck of the di-nuclear 
system is well formed. 

The fusion cross sections of ^^0+^''^Pb are also calculated with the ImQMD model by 
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FIG. 2: (Color online) Fusion excitation function of ^^0+^'^'^Pb reaction. The solid circles denote 
the experimental data [29]. The open circles denote the results of ImQMD with parameter sets 
IQ3 and SKP*. 

adopting IQ3 and SKP*. Through creating certain bombarding events (about 100) at each 
incident energy Ec.m. and at each impact parameter b, and counting the number of fusion 
events, we obtain the fusion probabihty gfusiEc.m.^b) of the reaction, by which the fusion 
cross section can be calculated [16]: 



o-fusiEc.m.) = 27r / bgfus db ~ 27r ^ fefi-fus A6. 



(9) 



The initial distance between the projectile and target is taken to be i? = 40 fm for calculating 
the fusion cross sections. Fig. 2 shows the comparison of our calculated results and the 
experimental data for the fusion reaction ^^0+^'^^Pb. The solid and open circles denote the 
experimental data and the calculation results, respectively. The measured fusion excitation 
function for ^^0+^°*Pb can be reasonably well reproduced with the new parameter set IQ3 
and SKP* at energies near and above the Coulomb barrier. The over-prediction of the fusion 
cross sections at sub-barrier energies is due to the shell effect of doubly-magic nuclei that is 
not well described with this semi-classical model. 

For further testing the reliability of IQ3 and SKP*, we have also calculated the charge 
distributions of fragments in multi-fragmentation processes at intermediate energy heavy- 
ion collisions. In Fig. 3 we show the charge distribution of fragments by using the ImQMD 
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FIG. 3: (Color online) Charge distribution of fragments for central collision of ^^ Au+^^ Au and 
^'^Ca+^'^Ca at 35AMeV. The solid squares denote the experimental data taken from Refs.[30, 31]. 
The open circles and open squares denote the results of ImQMD with IQ3 and SKP*, respectively. 

model with parameters set of IQ3 and SKP* for ^°Ca+^°Ca [30] and i^^Au+^^t^^ [g^j ^^ 
incident energy of 35 MeV/nucleon. Here we create 500 events for head-on collisions and for 
each event we self-consistently simulate the whole collision process till t = 6000 fm/c with 
a step size of At = 1 fm/c. We have found that the experimental data can be reproduced 
remarkably well. 

C. Dynamical nucleus- nucleus potential 

Based on the dynamical densities of the reaction system, the nucleus-nucleus potential 
can be obtained with the ETF approximation for the kinetic energies [15]. After the di- 
nuclear system is formed, the nucleus-nucleus potential may be described by a way like the 




FIG. 4: (Color online) Dynamical nucleus-nucleus potential for the reaction ^^O+^'^^Pb with pa- 
rameter set IQ2, IQ3 and SKP* at incident energies of (a) -E'c.m. = 80 MeV and (b) E'c.m. = 120 
MeV, respectively. Here we set the initial distance between two nuclei as 30 fm. The solid circles, 
open circles and triangles denote the results with IQ2, IQ3 and SKP*, respectively. The dashed 
line is to guide the eyes. 
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FIG. 5: (Color online) The same as Fig. 4, but for the reaction ^'^Ca+^'^'^Pb at incident energies 
of (a) Ec.m. = 179 MeV and (b) Ec.m. = 205 MeV, respectively. 



entrance channel potential [10] 

V{R)=E,,,{R)-E^-E2, (10) 

where EtotiR) is the energy of the composite system which is strongly dependent on the 
dynamical density distribution of the system obtained with the ImQMD model, Ei and 
E2 are the time average of the energies of the projectile and target nuclei, respectively. 
In this work, the dynamical nucleus-nucleus potential are calculated the same as in Ref. 
[15], but with the parameter sets IQ3 and SKP*. The kinetic energy coefficients cq, ci and 
C2 listed in Table II are determined by fitting the obtained kinetic energies T of a series 
of nuclei from light to heavy nuclei at their ground state with Eq.(7). To investigate the 
dynamical nucleus-nucleus potential, the head-on collisions of ^^0+^°^Pb and ^^Ca+^''^Pb 
at two different incident energies with three different parameter sets have been studied. As 
mentioned previously, the values of the incompressibility coefficient are the same but the 
wave-packet widths are different for IQ2 and SKP*, while the incompressibility coefficients 
are different but the wave-packet widths are the same for IQ3 and SKP*. 

Fig. 4 and Fig. 5 show the dynamical nucleus-nucleus potential for ^^0+^'^^Pb and 
'^^Ca-|-^°'^Pb with the parameter sets IQ2, IQ3 and SKP*, respectively. From the figures we 
find that: (1) the dynamical barrier height depends on the incident energy as mentioned 
in Ref. [15]; (2) the wave-packet width infiuences nuclear surface diffuseness and thus in- 
fluences both the potential barrier height and the potentials at short distances; and (3) 
the nuclear matter incompressibility seems just to affect the potentials at short distances if 
taking the same wave-packet width. Comparing the results with SKP* and IQ3, one sees 
that the potentials at short distances increase with the increase of the incompressibility co- 
efficient. To illustrate this point, we also study the static entrance channel potential [10] of 
^^Ca+^'^'^Pb with the Skyrme energy-density function by adopting different parameter sets. 
Fig. 6 shows the nuclear potential (i.e., removing the Coulomb potential from the entrance 
channel nucleus-nucleus potential) as a function of distance between two nuclei. Obviously, 
the nuclear potentials do increase with the value of K^o- In addition, the obtained dynamical 
nucleus-nucleus potential has been checked by directly using the barrier penetration calcu- 
lations for the fusion cross section. As an example, the fusion cross section for the reaction 
i6Q_|_208p|^ at Ec.m. = 80 MeV is calculated based on the obtained potential for head-on col- 
lisions with IQ3 which is shown in Fig.4(a). We find that the obtained fusion cross section 
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FIG. 6: (Color online) Nuclear potential of ^^Ca+^'^^Pb with different Skyrnie forces. The solid 
squares, stars and open circles denote the results with SKP, SkM* and SLy7, respectively. 

with the barrier penetration approach is very close to the result with Eq.(9). We also note 
that the obtained distribution function for the fusion probability g{ns{b) is different with the 
two approaches for this reaction, although the calculated fusion cross section is close to each 
other. The difference is due to that the reduced mass /i in the traditional barrier penetration 
calculations is fixed, but it changes as a function of distance R between two nuclei in the 
ImQMD simulations when the neck of the composite system is formed. 

To further investigate the influence of nuclear repulsion on the nucleus-nucleus potential, 
we study the nucleus-nucleus potential with IQ3 but varying the parameter /3. The value 
of P and the corresponding kinetic energy coefficients and Koo are listed in Table III (here 
C3 = 0). Fig. 7 shows the dynamical nucleus- nucleus potential and the contribution of the 
corresponding effective interaction potential energy and that of the corresponding kinetics 
energy for the reaction ^^O+^'^^Pb with different values of /3. One can see from Fig. 7 
(b) that the contribution of the effective interaction potential energy increases with the 
nuclear incompressibility coefficient as we expected. But we also note that the change of 
/3 causes larger change of the corresponding kinetic energy than that of the interaction 
potential energy, which results in the decrease of the total nucleus-nucleus potential at short 
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FIG. 7: (Color online) (a) Dynamical nucleus-nucleus potential, (b) contribution of the interaction 
potential energy and (c) contribution of the kinetics energy for reaction ^^0+^'^'^Pb at -Ec.m. = 80 
MeV based on the parameter set IQ3 but varying the value of /3. The solid circles, crosses, open 
circles and triangles denote the results with /3=138, 140, 142 and 145 MeV, respectively. 

TABLE III: The same as Table II, but taking different value for /3. 



Parameter 


Co 


(MeV /m2) 


Cl 


(MeV /m2) 


Koo (MeV) 


/3 = 138 




43.3 




0.34 




225 


/3 = 140 




40.8 




0.30 




229 


/3 = 142 




38.0 




0.83 




232 


/? = 145 




33.8 




1.43 




237 



distances with increasing the nuclear repulsion. 

IV. CONCLUSION 

In this work, the dynamical nucleus-nucleus potentials for fusion reactions have been 
investigated by using the improved quantum molecular dynamics model with different sets 
of parameters. By using two new sets of parameters IQ3 and SKP* with which the measured 
fusion excitation function for ^^0-|-^°^Pb and the charge distribution of fragments for Ca-|-Ca 
and Au-|-Au in multi-fragmentation process can be reasonably well reproduced, we find that 
both the nuclear incompressibility and the range of nucleon-nucleon interactions significantly 
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influence the nucleus-nucleus potential. The interaction range represented by the gaussian 
wave-packet width affects both the barrier height and the potentials at short distances. 
The incompressibility coefficient of nuclear matter mainly influences the potentials at short 
distances if taking the same wave-packet width. In addition, the nuclear repulsion influences 
both the effective interaction potential energy and the kinetic energy of a fusion system. 
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